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BOUNDARY DRIVEN REACTION DIFFUSION MODEL 
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Abstract. We derive a large deviation principle for the empirical currents of lattice 
gas dynamics which combine a fast stirring mechanism (Symmetric Simple Exclusion 
Process) and creation/annihilation mechanisms (Glauber dynamics). Previous results on 
the density large deviations can be recovered from this general large deviation principle. 
The contribution of external driving forces due to reservoirs at the boundary of the system 
is also taken into account. 

1. Introduction 

A major challenge in non-equilibrium statistical physics is to provide a counterpart to the 
notion of free energy and to the Gibbs measure which is the cornerstone of the equilibrium 
statistical physics. Large deviation principles have been proposed as a good alternative 
to investigate properties of non equilibrium systems [G] 01 BDGJL3]. In particular, a 
lot of attention has been devoted to the case of lattice gas dynamics for which explicit 
large deviation principles can be derived in the hydrodynamic scaling (we refer the reader 
to [Dj [BDGJL3J for recent surveys and further references). Motivated by these recent 
progress in physics, the original mathematical works |KOV[ IDV] on the hydrodynamic 
large deviations for conservative dynamics have been generalized to take into account the 
contribution of reservoirs at the boundary of the system [BDGJLTJ IBLM[ IFLMj and the 
current of particles flowing in the system [BDGJL2J. 

An interesting class of models has been introduced in |DFLj to describe reaction diffu- 
sion equations by combining the Symmetric Simple Exclusion Process (SSEP) to a Glauber 
dynamics which models the annihilation and creation of particles. In |DFL] . the hydro- 
dynamic limit as well as the fluctuations of the density have been investigated for these 
models. The density hydrodynamic large deviations have then been proved in (JLV] . In 
this paper, we generalize this result by deriving the joint large deviations of the density 
and of the (conservative and non-conservative) currents flowing in the system. We also 
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take into account the contribution of reservoirs acting at the boundary of the system. 
Our results were motivated by the recent research in non-equilibrium statistical physics on 
dissipative dynamics and in particular on granular media [El ILL] ISL] . We refer to |BLj 
for a more comprehensive discussion on the physical aspects of the large deviations for 
dissipative systems. 

Contrary to the purely conservative dynamics [BDGJL2J, one has to introduce two types 
of currents: the conservative integrated current Q t which records the particle jumps from 
the diffusive part of the dynamics (SSEP) and the non-conservative integrated current K t 
associated to the creation annihilation process (Glauber). Heuristically, if one denotes by 
Qt( r ) an d K t (r) the instantaneous currents at time t and location r, then the density obeys 
the following equation 

d tPt (r) = -d r Q t (r) + K t (r) . 

The goal of this paper is to compute the asymptotic cost of observing an atypical trajectory 
of the currents Q,K and of the density p when the number of particles tends to infinity. 
Even so it is apparently more complicated to consider the joint deviations of the empirical 
currents and the empirical density, it turns out that the structure of the joint large deviation 
functional Iq is more transparent as it splits into two distinct contributions involving either 
the diffusive part or the Glauber part of the dynamics 

I (p,Q,K) = h{p,Q)+I 2 (p,K). 

The precise form of the functional Iq can be found in (I2.12p . In section 0, the density 
large deviation functional derived in |JLVj is recovered by a contraction principle. This 
provides a natural interpretation of the density large deviation functional as the optimal 
combination between the two macroscopic currents Q, K in order to create the atypical 
density trajectory p at a minimal cost I\ + 

Our proof relies on the standard machinery developed to study hydrodynamic large de- 
viations |KOV[ IKLj . as well as on more recent tools introduced in BDG.ILU IBDGJL2"| 
IBLM[ IFLMj . Therefore in this paper, we will not detail the aspects of the proof which can 
be deduced readily from the existing literature and we will focus on the new features oc- 
curring from the non-conservative part of the dynamics. The paper is organized as follows. 
In section [21 we introduce the model and state the main results. A strong form of local 
equilibrium is stated in section [3j The upper and lower bound of the joint density/ current 
large deviations are derived in sections [4] and [5j Finally the density large deviations are 
recovered in section El 

2. Notations and Results 

2.1. The microscopic dynamics and Notations. We consider the one- dimensional 
Symmetric Simple Exclusion Process (SSEP) in the domain {— N, ...,N} with creation 
and annihilation of particles in the bulk and reservoirs at the boundaries. More precisely, 
the particles perform random walks with an exclusion constraint which imposes at most 
one particle per site and particles can be removed or created in the bulk according to a rate 
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which depends on the local configurations. At the boundaries ±iV, two reservoirs maintain 
constant densities. The dynamics can be viewed as a toy model for chemical reactions 
where the chemicals are injected at the boundaries, then diffuse and react in the system 
[DFL] . 

The stochastic dynamics is a Markov process on {0, 1} 27V+1 whose generator is obtained 
by adding the generators of the different dynamics 

N 2 N 2 N 2 

(2-1) L N = -^-L 0jN + -^-L +i jy + —^-L_ N + Li^ N , 

with the SSEP generator 

N-l 

wfo) = E [f(v x ' x+1 ) - m] , 

x=-N 

and creation and annihilation generators at the boundaries depending on the parameters 
P+ and p~ 

L+, N f(v) = [V(N) + - rj(N)] [f( v N ) - f( v )] 

L-, N f(v) = H-N) + - v(-N)} [f( V - N ) - M] 

where 

/ \ -r , { v(y) if z = x 

r)(z) it z x xv/ \ ) ) \ -c 

J f \ . £ ; ri x ' y (z) = < Ti(x) if z = y . 

1 - T)(z) if Z = X ' ' U I / \ i 

^ 77(2) else 

Finally the creation and annihilation generator in the bulk is given by 

AT-Af-l 
x=-V+M+l 

where the rate of creation and annihilation c(x, ■) is a non negative cylindric function with 
range M, i.e. there exists a fixed integer M (possibly equal to 0) such that c(x, rj) = 
c(Vx-m, ■ ■ ■ 1 Vx+m) depends only of the values of rj in {x — M, . . . , x + M}. Remark that, 
the diffusive part of the process is speeded up by iV 2 to obtain a non trivial hydrodynamic 
evolution. 

For a given trajectory rj : [0, T] — > {0, 1} 2N+1 , let pf be the empirical density of particles 
in [-1,1] at time t £ [0, T] 

1 N 

Pt = n E Vt(x)5 x /N- 

x=-N 

We denote by (x) the conservative current through the edge (x,x + 1), i.e. the total 
number of particles that have jumped from x to x + 1 minus the total number of particles 
that have jumped from x + 1 to x between the times and t. The empirical measure 
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associated to this current is denned as the signed measure on [—1, 1] 

N-l 

(2-2) Qf = ^ Yl Q?(*)8*/N' 

x=-N 

The renormalization by N 2 takes into account the space renormalization as well as the 
diffusive scaling of the SSEP dynamics which leads to an extra factor N. We denote by 
K^(x) the non-conservative current at site x, i.e. the total number of particles created 
minus the total number of particles annihilated at site x between times and t. The 
corresponding empirical measure is 

1 N 

x=-N 

For any continuous function if G C([— 1, 1]), we will use the notation 

JV-l 
x=-N 

The same notation will be used for . As the conservative current applies to edges, we 
will write for ip G 1, 1]) 

x=-N \ \ / / 

Note that <p (^jf) —f is of order 1/N so that the scaling is coherent with (12. 2p . Finally, 
for any functions f(s,r) in [0, T] x [—1, 1], we use the shorthand notation 

VsG[0,T], (/„) = ^ drf(s,r). 

The density profiles bounded away from and 1 will be relevant so that we introduce 
C e ([— 1, 1]) the set of continuous functions / on [—1,1] for which there exists a constant 
e > such that e < / < 1 — e. Given a function 7 G C e ([— 1, 1]), let u? be the Bernoulli 
product measure on {— N, N} with marginals 

(2.3) VfcG{-iV,iV}, <fo(*) = l)=7(^). 

For a G [0, 1], the Bernoulli product measure with uniform density a is denoted by u^. Let 
be the probability measure associated to the Markov process (p^ ,Qf , K^) on [0, T] 
with initial measure on the particle configurations. 

We define M. the set of signed measures on [—1, 1] endowed with the weak topology. We 
also consider Ai the subset of M. of all absolutely continuous measures wrt the Lebesgue 
measure with positive density bounded by 1 

M = {p(x)dx G M, < p(x) < 1 a.e.} . 
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In order to consider the joint large deviations of (pf ,Qf , K^) during the time interval 
[0, T], we will work on £ = D{[0, T], Aio xMx M) the space of cad-lag trajectories with 
values in A4q xM xM endowed with the Skorohod topology [EKJ. 

2.2. The results. The hydrodynamic behavior of the microscopic dynamics introduced 
in section I2.1[ can be described in terms of a few macroscopic parameters 

(2.4) Va G [0, 1], C(a) = u a (c(0,r])(l - 7/(0)), A(a) = ^(0(0,77)77(0)) , 
where C and A represent the average creation and annihilation rates at density a and 

(2.5) VaG[0,l], (7(a) = u a ( Vo (l - m )) = a(l - a) , 



which is the conductivity of the SSEP. Finally, we denote by p± 
the boundaries imposed by the reservoirs. 

We are now ready to state the hydrodynamic limit 



1 + & 



the densities at 



Theorem 2.1. Let 7 G C e ([— 1, 1]) with 7(— 1) = p~ and 7(1) = p + . For each T > and 



V<5 > 0, 



lim P 

JV-s-oo 



A' 



dt(pt ft) ~ {pt ft) 



> 5 



0. 



where p(t, x) is the unique weak solution of 

dtp(t,x) = ±Ap(t,x) + C(p(t,x))-A(p(t,x)) 



(2.6) 



p(t,±l) 
p(0,x) 



P± 
7(x) 



The meaning of weak solution of (12.61) is recalled in the Appendix (with H = G = 0). 
As a consequence of Theorem 12.11 a law of large numbers holds for the currents. 
Theorem 2.2. For any test function (p G C 1 ([— 1, 1]) and t > 

1 



V5 > 0, 



lim P. 

Af-s>oo 



N 



(Q?p) + 



lim P. 



N 



(K t N <p)~ 



ds (<f Vp s ) 

ds (f (c(p s ; 



> 5 







A(Ps)) > 



> 5 



where p is the unique weak solution of $2.6}) . 



Theorems 12. 1[ 12.21 can be deduced from the methods used to prove the large deviations 
so that their derivation is omitted (see also |DFL[ IBDGJL2"] ). 

Before stating the large deviation principle, we need more notation. Let G and H be 
smooth functions in [0, T] x [—1, 1]. For a given trajectory (p, Q, K) in £, we set 

(2.7) J G , H (p, Q, K) = 4(p, Q) + (p, K) 



THIERRY BODINEAU AND MAXIME LAGOUGE 



with 

(2 ' 8) 

J 1 H {p,Q)={Q T VH T )- J (Q a JLvH a ^ds-~^ ( Ps AH s )ds 

< a ^)( V ^) 2 > ds + \p+ [ X ) ds ~ \p- [ - 1 ) ds > 

where o~{u) = u(l — u) is defined in (12.51) and 

(2 ' 9) 

J 2 G (p,K) = (K T G T )- J {K a ^-Gs\ds- J ds (C(p s )(e G °-l) + A( Ps )(e- G °-l)) , 

where A and C were introduced in (12.41) . 

The first functional is related to the contribution of the conservative currents and the 
second one to the non-conservative currents (see theorem 12.41) . We define 

J (p, Q, K) = sup Jg,h (p, Q, K) = sup Jh(p,Q)+ sup J%(p,K), 

G,H HeC 1 ' 2 GeC 1 ' 

where the supremum is taken on regular functions G and H . Note that the functions G 
and H can take arbitrary (finite) values at the boundaries. 

Define A as the set of trajectories (p, Q, K) satisfying the following two conditions: 

• Conservation law. For all test function (p G C 1 ([— 1, 1]) vanishing at the boundaries 

(2.10) (p t ip) - (p oV ) = (Q t d x <p) + (K t <p), Q = 0, K = 0. 

• Energy condition. The energy Q(p) of the density trajectory is finite with 



Q(p) = sup | / dt j dx p(t, x)V(p(t, x) / dt I dx 

f kJo J-i 2 J J_ l 



<p(t,x) 



(2.11) - f dt(p +V (t,l)-p^(t,-l)) 

Jo 

and the supremum is taken over smooth functions (p in [0,T] x [—1, 1]. 
For smooth trajectories the conservation law reduces to 

dtp + d x d t Q - d t K = , 

where d t Q, d t K are the instantaneous currents and the energy condition reads 



Q (ti = lJ o dt J dx (Vp(t,x)) 2 <oo. 



The energy condition was introduced in |QRV[ IBLMt IFLM] to control the approximation 
procedure in the derivation of the large deviation lower bound (see Theorem 12. 6p . 
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Finally we define the dynamical rate function 
(2.12) IMQ,K)-{ J }^ 1^^' 

To take into account the large deviations of the initial measure , we introduce for any 
function m : [— 1, 1] — >■ [0, 1] 

/ 771 \ I 1 — 772 

h 7 (m) = ( mlog — ) + ( (1 — m) log 



7 / \ 1-7 

and 

h^im) if p,(dx) = m(x) dx , 



j , x _ f h^(m) if p(dx) = 
7 ^ y +oo otherwise . 



The rate function is the sum of the dynamical deviation cost from the initial measure 
and the deviation cost from the hydrodynamic trajectory 

I(p,Q,K) = I (p,Q,K) + I 1 ( Po ). 

From now, the initial density profile 7 is a given smooth function in C e ([— 1, 1]) equal to 
p± at the boundaries. We state the large deviation Theorems. 

Theorem 2.3. For all closed set T e S 

limsupllogP? [(p»,Q?,K») 6 7] < - inf I(p,Q,K). 



We first state the lower bound for regular trajectories. 

Theorem 2.4. Lei (p, Q, if) be a regular trajectory. Then the large deviation functional 
I has an explicit form 

f T f 1 ( (d t Q(x,t) + d x p(x,t)) 2 / , x 

(2.13) I (p,Q,K) = j( dt j( dx j 2a{p(x t)) + ^{p(x,t),d t K(x,t)j 

with 

(2.14) $(p, «) = C(p) + A(p) - V« 2 + 4A(p)C(p) + .log ( ^ 2 + yj^ (P)+ ~) 

IfC(p) = 0, then <£> becomes 

'^(p) + «-«log(^), if «<0, 



$(p, «) 



00, if k>0. 



^4 similar formula holds if A(p) = 0. 

For any open set O G £ containing the regular trajectory (p, Q, if) 

liminf llogP^ [(pf ,Qf ,iff) £0]) -I(p,Q,K). 



8 



THIERRY BODINEAU AND MAXIME LAGOUGE 



Remark 2.5. The first contribution to ( 12 . 1 3[) comes from the difference between the in- 
stantaneous empirical current d t Q and the canonical instantaneous current associated to 
p which is — |Vp. This term has already been analyzed in the conservative dynamics 
[BP I IBDGJL2"] . The second term in (12 . 131) should be interpreted as the large deviation 
functional associated to Poisson processes with parameters C(p t ) and A(p t ). 

In order to derive the lower bound for general trajectories, we introduce two technical 
assumptions on the rates (12. 4ft : 

Assumption (LI). The rate A (resp C) is either concave and positive on ]0, 1[ or uni- 
formly equal to zero. 

Assumption (L2). The functions A and C are monotonous and 
(2.15) G [0, 1], A'(z) > 0, C\z) < . 

Then 

Theorem 2.6. Assume (LI) and (L2), then for all open set O G S, 

hminf llogP^ [(p»,Q?,K?) 60]>- inf I(p,Q,K). 

N-+oo iV (p,Q,K)eO 



The concavity assumption (LI) has been introduced in [JLVJ. As we shall see in Section 
15.21 it simplifies the proof of the lower bound, however it is mainly technical and Theorem 
12.61 should be valid without assumption (LI). We refer to |QRV, BLMl IFLM] for further 
results on this generalization in the case of conservative dynamics. Assumption (L2) will 
be used in the Appendix only to ensure the uniqueness of the weak solutions for singular 
perturbations of the hydrodynamic equation (12.61) . 



3. Modified dynamics and local equilibrium 

Local equilibrium lies at the heart of the hydrodynamic limit theory and it states that 
during the time evolution the local measure remains close to an equilibrium measure with 
a varying density. In this section, we state, in our framework, a strong form of local 
equilibrium which will be useful for the derivation of the hydrodynamic large deviations. 
The proofs are omitted as they follow the scheme introduced in [KOVl IKL| IBDGJLl] . 



3.1. The modified dynamics. We first define a modification of the original process (12. ip 
which will be used to derive the large deviations. For regular functions G and H on 
[0,T] x [—1, 1], denote L c ^ )Ht the time dependent generator given by 

fqi\ T Gt,H t _ N 2 H Q t N 2 , N2 r 

(3.1 J N - -^ h 0,N + L l,N+ ~^ L +,N + ~2 L -< N ' 
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with 



N-l 



E if(v x ' x+1 ) - f(v)] exp [( V (x) - V (x + 1)) (H{t, ^±i) - , j-] 

=-N 
N 

E c (^) h(^)e" G(WAr) + (1 -^(x))e G ^] [/(,/■) - /fo)] . 



x=-7V 



The modified dynamics induces a weak drift in the conservative dynamics. For large N, 

a particle jumps from x to x ± 1 at rate — ^1 ± —d x H (t, ff)J ■ L Gt N is the generator of 

a non-conservative dynamics for which the intensity of creation and annihilation varies in 
time and space according to G. Finally, let W^ t G,H ^ e the probability measure associated 
to the process with initial measure and generator L N U *. We stress the fact that the 
reservoir dynamics are unchanged. As in Theorem I2.1[ one can show that the modified 
dynamics follows the hydrodynamic limit equation 

(3.2) d t p(t,x) = ±Ap(t,x)-V(a(p(t,x)) VH(t,x))+C(p(t,x))e G ^-A(p(t,x))e- G ^ . 



3.2. Local equilibrium. We set 

define the local density as 



{y G {— N, . . . , N} such that \y — x\ ^ 1} and 



(3.3) 



\ A i 



yeA 



(x) 



where \A\ | stands for the number of sites in A\ x> . Let ip be a cylindric function with 
support in {—R, . . . ,R}. Given ip a regular function on [0, T] x [—1, 1] and 5, e > 0, we 
define the set B$ etp (ip) on the trajectories {t} s } s <t as 
(3.4) 

„T N-R 



(x) 



Bs,e,<p{lp) = < {Vs}s<T, 



ds 



1 



2(iV - R) 



E ^ ( s ' a?) [^( r ^) - Vc*)^)] 



x=-N+R 

where T x rj is the configuration rj shifted by x. 

The local equilibrium property also holds for the modified dynamics (13. ip . 

Theorem 3.1. Given cp, ip and 6 > 0, the trajectories concentrate super- exponentially fast 
on the set Bg >S)tp (ip) 



lim sup lim sup — logP G H [ Bs^W 



-OO . 



Moreover, the reservoirs impose local equilibrium at the boundaries with the densities 
p + ,P-. For any continuous function <£> in [0,T] 

-T 



lim sup — logP G H 



ds$(s)( Vs (±N)-p ± ) 



> 5 



— oo. 
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The derivation of Theorem 13.11 follows from the bounds on entropy production |KOVj 
which can be adapted to control the boundary terms as in [BD GJLi] . 



4. Large deviation upper bound 

The derivation of Theorem 12.31 is split into several steps. First an upper bound with 
the rate function J (I2.12p is derived for compact sets, then for closed sets. Finally, we 
prove that the rate function is infinite for the trajectories which do not belong to the set 
A introduced in ( l2TT0j) . (IXTTD . 

4.1. The upper bound for compact sets. In order to compare the original dynamics 
starting from the initial profile 7 to a modified dynamics with regular drifts G, H ( 13. ip 
starting from the initial profile u, we compute the Radon- Nikodym derivative |KL] 



,G,H 



exp 



n{(Q$VH t )-J ds(Q? ±VH S ) + (K»G T )- J ds(K?±G s ) 
(4.1) +0 H (^))- j\ N (s)ds + h^(Po) 



where h 7i0J (p) = (p \og(u/j)) + ((1 - p) log((l - w)/(l - 7))) and 



N 



N-l 



Rn(s) = — Yl Vs(x)(l- Va (x + 1)) 



x=-N 



exp [His 



x + 1 
N 



+7] s (x + 1)(1 - rj s (x)) 



exp ( H (s, 



H s 



x + 1 
N 



N 



(4.2) + c ( x ^s) [exp ((1 - 27?.0r)) G(s, -)) - 1 

x=-N 

expanding the exponential and summing by parts, we get 

2V-1 r 



Rn(s) 



1 E Vs(x)AH(s^) 

[_ X =-N L 

+- (VH (s, ^ 2 [ Vs (x)(l - Vs (x + 1)) + Vs (x + 1)(1 - Vs (x))} 



+-\ Vs (N)VH(s, 1) - rj s (-N)VH(s, -1) 



N 



+ c ^ Vs) [exp ((1 - 2 Vs (x)) G(s, -|)) - lj + 0(1) . 



x=-N 
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Thanks to the local equilibrium, the microsopic expressions in Rn(s) can be replaced by 
their averages. We set 



Bs,e = £«, e ,(Vfl) a (77o(l - rjo)) f]B 5 ,e,(cG-i)(c(0,r])r] ) P|5f ie) ( e -o_i)(c(0,7;)(l - rj )) f)B' s 
(4.3) 

where the sets Bg £jV3 (-) were introduced in (13.41) and 



± ■ 



B' 



5,± 



[ dsVH{s,±l)(ri 8 (±N)-p±) 
Jo 



^ 5 



The super-exponential replacement Theorem 13.11 implies 



V<5 > 0, limsuplimsup — logK (Bg e ) = —00. 
For trajectories in Bs )£ , the Radon-Nikodym derivative can be approximated as follows 



dP N 



G,H 



cxp 



N (j G , H ,e{p N , Q N , K N ) + h 7> M + O g>h (±) + 0(S^j 



where we used the functional 
(4.4) 

J GiH!e (p,Q,K) = (Q T VH T )-fids(Q s £ s VH s ) + (K T G T )-fids(K s £G s ) 



1 r T 



2 J ds(p s AH 8 + <7(p 8 *i £ )(VH s y) 
- J ds (C(p s * i e )(e G ° - 1) + A(p s * i e )( ( 



-G 3 



i)> 



+i / 5 + / T <isVff( S , 1) - ip_ J*dsVH(s, -1) , 
with A, C and cr as in (12 .4p . (12. 5p . The function i t is an approximation of unity 



Ps * i e — f/ s . 



For any (9 open set in £, we can then deduce the large deviation upper bound 



limsup — P^(£>) < max 



limsup - logP^ ((p N , Q N , K N ) e O H B s , £ ) ; 

N^oo JV 



limsup- log py(B^ 



^ limsup — log W% GH 

N^oo JV 



< limsup — log P^ G)i y 



dP 



dP N 

" 7 1 fjv zv^ 



iV->oo N 



-N(j G ,HA(p N ,Q N >K N ))+h-,, ul (p^)+OG,H{jr)+°( s )) l 



to 



SUp I - J G ,H,e{P> Q> K ) ~ KAPo) }+0(S). 



(p,Q,K)eO 
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This is true for any (e, G, H, u) and any 5 > so finally 
1 



(4.5) limsup^-P^(C) ^ inf 



sup 



JG,H,e(P: Q, K) - h JjU (p ) >. 



The previous bound can then be extended to any compact set K by using a finite covering 
with open sets (see [KLJ) 



limsup ^-logP^(/C) < - sup sup{J Gi #(p, Q, K)} = -J(p,Q,K). 



N 



{p,Q,K)elC G,H 



4.2. The upper bound for closed sets. We are going to prove the exponential tightness, 
i.e. to exhibit a sequence {JC n } of compact sets in £ such that for any n 



(4.6) 



limsup-logP; v (/q) < -n. 



iV-i-oo 



The large deviation upper bound will then follow for general closed sets F of £ by noticing 
lim sup 1 log P^ (F) ^ lim sup 1 log [F% (F n K n ) + P?(/C°)] 



< max 



sup J(p, Q, K) ; — n 

(p,Q,K)£Fn)C n 



Letting n go to oo completes the upper bound for closed sets. 

In order to build a sequence of compact sets, we need to check first that the measures 
concentrate on equi- continuous trajectories 

Lemma 4.1. Let <p be a C 1 function on [—1, 1]. Then we have for all e > 0, 



(4.7) 
(4.8) 
(4.9) 



lim lim sup — log P, 

<5->0 at^oo N 



sup I <pf 0> - <pf 0> I > 

|t-s|<<$ 



lim lim sup — log P 

5^0 tv^oo N 



lim lim sup — log P, 



<5->0 TV 



N 



7 



sup I (Q? 0) - (Q? cf>) | > e 

|l-s|<<$ 

sup | (KU) - (K? <f>) | > e 

|t-s|<5 



— OO 



— OO 



-oo . 



where the supremum is taken over s,t in [0,T]. 

Second, we need estimates on the total variation norm of the empirical currents 
Lemma 4.2. For any time T > 0, one has 



(4.10) 



lim lim ^logP^ (sup 1 \ K ?( X )\ ^ a ) 

a-+oc N^oo Jy ' \ o<i<T -/V z — ' 1 1 / 

\ x=-N / 



-OO . 
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(4.11) 



N 



lim lim — log P^ 

a-*oo AT->oo N 1 



f sup 1^2 S \Qt(*)\ > a ) 

\0<t<T iV 2 x ^ N ) 



-oo . 



Should the currents be bounded (as the density), then Lemma [4.11 would be enough to 
ensure the exponential tightness |KL] . 

We postpone the derivation of the Lemmas and conclude the proof of the exponential 
tightness (14.61) . We focus on the conservative current as the same strategy applies to p, K. 
Fix a sequence (pi of C 2 functions on [—1, 1] dense in C([— 1, 1]). We set 



Q G D([0,T],M) 



sup 

|*— s|<<5 



(Qt<f>i)-(Qs<f>i) 



re- 



using Lemma [4.11 we have, 

V/ ^ 0, Vm ^ l,Vn > l,35(n,m,/), F%(Q N $ Ci An , mAm ) < exp(-Nnml) . 

We introduce also C' n = {Q G D([0,T],JA) ; sup is ; T \Q t \ ^ k(n)} , where \Q t \ stands for 
the total variation norm of the measure Qt and k(n) is chosen, according to Lemma [4.21 
such that 



Now consider 
(4.12) 



F»(Q N i O < exp(-iVn) 



fcn — f} Cl,6{n,m,l),m C' n 



Ascoli Theorem (see |EK] Theorem 6.3 page 123) implies that K n is a compact set. Com- 
bining the previous estimates we see that K n satisfies (14. 6p . 



Proof of Lemma \4-l\ We start by proving (14.81) and follow the strategy of [BDGJL2] . It is 
enough to show that the expression below goes to — oo as 5 vanishes 



max lim sup — log P. 



sup (Qf 0>-<Qf 5 0>>€ 

fc<5<t<(fc+l)<5 



Thanks to the computation of the Radon-Nikodym derivative of the modified dynamics 
(j42D (with VH — jjtfi, G — 0), we know from (@~T} that for all a > 0, 



{ i r T JV_1 / 

explN(Q% a( p)-- / ds ]T U(x)aV4>(£ 

I ^° x=-N \ 



1 

+ 4 







a<f> 


m 







[r) s (x)(l - rj s (x + 1)) + r} s (x + 1)(1 - rj s (x))) 
+N[r) s (N)<f>(l) - r] 8 (-N)<l>(-ij\ + 0(1) 
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is a mean one positive martingale. One easily checks that the integral term is bounded 
above by C^a 2 NT (we will take the limit as a goes to infinity), where is a constant 
depending only on 0. Therefore, multiplying by aN, adding and substracting the integral 
term of the logarithm of the martingale and exponentiating, 
(4.13) 



sup (Qf 0>-(Q^0)> e 



sup * > exp {eaN — C^a 2 N5} 

k6^t^(k+l)S Jvt-Sk 



Taking 5 small enough such that C^,a 2 N5 < ^p, then Doob's inequality implies that the 
last expression is bounded above by exp[— aiVe/4]. Letting a go to oo completes the proof. 
For the non-conservative current, (14. 9 j) follows in the same way by using the martingale 



Aft = exp 



(p4>K$) - / ds V c(x, Vs ) [ e ^-^))H-/N) _ X ] 

Jo „__ A r 



x=-N 

The proof of (14. 7p can be found in [KLj . □ 



Proof of Lemma \4-S\ The probability of the event in (14. lOj) can be estimated from above 



by the large deviations of (2iV + 1) independent Poisson processes thus 

(4.14) Va > 0, ( sup i- V \K?(x)\ > a ) < exp(-NC a ) 

\o<t<T N ^ N J 

where C a goes to infinity as a diverges. 

We turn now to the bound (14. lip and show that 

(4.15) Va > 0, P* ( sup -L \Q?(x)\ > a ) < e 



NCL 



0<t<T N 2 

x=—N 



where C' a goes to infinity as a diverges. To prove (14.151) . we first use a microscopic identity 
(which holds at any time) 

(4.16) V* E {-N, N}, Q» (x - 1) - Q»{x) = t£{x) - rff{x) - K t N (x) . 

Therefore (I4.14p implies that with probability at least 1 — e~ NCa , the conservative current 
through the edge (x, x + 1) satisfies 

Vt G [0,T], |Qf (x)\ < \Q?(-N)\ + 2N(1 + a) . 

From this, we get 

1 \| \ K„ 1 / „-NC a , tb>N ( \Q?{-*0\ 



7 



SU P m E l^f (*)| > 5a U e"^ + P W sup > 3a 

. o<t<T N 1 J~* N J \ 0<t<T N 



<^ e -Nc a t + F N( Q[(N) >3a _2] +f n( inf Qf( N ) < _ 3a + 2 ] 

7 \0<t<T N J 1 \0<t<T N J 

(4.17) 
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We bound now the first term on the RHS (the second one can be bounded similarly by 
symmetry). Using again the identity (I4.16p . we get 

AT-l 



P^fsup ^-^>3a- 2 )<pW sup -L £ Q»(x)>a-A]+e 

\0<t<T N J \0<t<T N 2 x t^ N I 



-NC a t 



The terms in the RHS can be estimated as in (14.131) by reducing to a martingale estimate. 
This completes (14.151) . □ 

4.3. The set A. To complete the derivation of the upper bound, we prove that the tra- 
jectories concentrate exponentially fast on the set A introduced in (I2.10p . (12.111) . 

The conservation law. 
Let Cq([— 1, 1]) be the set of twice different iable functions vanishing at the boundary. 

Lemma 4.3. For any <p G Cg([— 1, 1]), we introduce 

N-l 



x=-N+l 

Then for any 5 > 0, 



limsup^logP^ [\V T (p N ,Q N ,K N ,<p)\ >5}= -oo. 

Proof. For any site x and times s < t, the following microscopic relation holds 

n?(x) - r£{x) = [Cff(x - 1) - €ff(x - 1)] - [Q?(x) - (ff (x)) + K»{x) - K?{x) . 
Summing in x and integrating by parts the Q-term gives 



1 ^ r(x+l)/N , 1 

V t { P n ,Q n ,K n ,4>) = Tt Y j W x Q%{x), with W x = du( X 



x/N N 



uU"(u) 



[X 



x=-N 

where we used that <p vanishes at the boundaries. From this identity we get for any a > 
(4.18) P^ {V T (p N ,Q N ,K N ,^ >6)^ e- aNS P^ (exp (a £ W X Q N T ' 

\ \ X =-N 

We first note that W x is of the order 1/N 2 and that uniformly in x, W x+ \ — W x is of the 
order e^/N 2 where vanishes to as iV goes to infinity. Thus the identity (14.11) implies 

oT N ( N ~ x \ 

= exp [a W x Q»(x) + ao^N) \ , 

7 \ x=-N J 

where H is such that H ((x+l) / N) — H (x / N) = W x . This function is a mean one martingale 
and (I4.18P leads to 

1 \ogF»(v T (p N , Q N , K N , 0) > 5) ^ -aS + ao^ (1) . 
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Letting iV and then a go to oo concludes the proof of the lemma. □ 

Let (<pk) be a dense sequence of functions in C5 3 ([— 1, 1]) and define 

A s n = {(p,Q,K) G £ such that Vfc<«, \V T (j>,Q,K,<j> k )\ < 5}. 
The previous lemma gives that for all n and 5 > 0, 

limsup ^ logP^ [(p N , Q N , K N ) GF]<- inf J(p, Q, X). 

This is true for any n and (<j) n ) is dense. So letting 5 go to 0, we can define Iq = +oo for 
the trajectories which do not satisfy the relation (I2.10p . 

Energy condition. 

Lemma 4.4. For any smooth ip : (0, T) x [—1, 1] — > M.. There is a constant cq such that 

(4.19) VA; > 0, lim sup limsup — log ( Q v (p N ) > k) < -k + c T , 
where Q v is defined for some suitable constant c 

QM= [ dt f dx p(t,x)Vcp(t,x) - [ dt(p+<p(t,l) -p_p(t,-l)) 

(4.20) - — I dt I dx<p{t,x) 2 . 

2c Jo 7-i 

The proof of this Lemma follows from |BLMl IFLM] and therefore is omitted. By consid- 
ering a dense sequence of functions {y?fc}, one can deduce from (10511 [FLM] that the large 
deviations are infinite for the trajectories p such that sup^ Q<p(p) = +oo. Note that if p is 
such that sup^, Q v (p) < +oo, then p is in L 2 ([0, T], H(] — 1, 1[)) and Riesz representation 
theorem implies that there exists Vp such that for any (p : (0,T) x [—1, 1] — > R. 

/ cit / dxp(t, x)V(p(t, x) — {p + (p(t, 1) — p~-<p(t, — 1)} = dt dxV p(t, x)tp(t, x) . 
Jo J-i Jo J-i 

(4.21) 

5. Large deviations lower bound 

The derivation of the lower bound is split into two parts following the general scheme 
for hydrodynamic large deviations [KOVf |KL] . First we derive the lower bound for regular 
trajectories. Then under assumptions (L1,L2), we prove that general trajectories can be 
approximated by regular trajectories. A key feature in this approximation procedure is 
that the contribution of both currents decouple 

(5.1) I q (j>,Q,K)=I 1 (j>,Q)+I 2 (j>,K)= sup J 1 H (P,Q)+ sup J&p,K), 

where Jjj and Jq were defined in (12.81) and (12. 9p . This simplifies some steps in the ap- 
proximation as both functionals can be analyzed independently contrary to the study in 
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[JLV] . Assumption (LI) provides some convexity properties of the large deviation func- 
tional which simplify the proof. We follow closely some arguments of |KL[ IBDGJL2] . Thus 
we will sketch the main steps of the proofs and only detail the new aspects related to the 
non-conservative currents. 



5.1. Lower bound for regular trajectories. In this section we derive Theorem 12.41 
Suppose p, Q, K are regular in time and space and that for any time t and x, p(t, x) is 
bounded away from and 1. To the trajectory (p,Q,K), one can associate the functions 
G(t,x) and H(t,x) satisfying 



1, 

2 

d t K{t,x) = C(p(t,x))e G ^ - A(p(t,x))e- G ^ 



(52) i d t Q(t,x) = -~Vp(t,x)+<r(p(t,x))VH(t,x), 



The pointwise existence of G comes from the fact that the polynomial CX 2 — KX — A 
has one positive root for positive C and A. Furthermore H is well defined as long as 
a(p) = p(l — p) 7^ 0. We choose H(t, —1) = and the value of H at x — 1 is imposed by 
equation (15. 2p . contrary to the case of density large deviations where H is equal to at 
both boundaries (see section [6]). 

Lemma 5.1. For regular trajectories (p,Q,K) the functional 1$ $2.12\) is given by 
(5.3) I (p,Q,K) = J G>H (p,Q,K), 

with G and H as in Ii5. || ). This expression coincides with the explicit form of the functional 

(EH. 

Proof. As (p, Q, K) is regular in time, Jq ^ can be rewritten after integration by parts as 



(5.4) J^ R (p,Q,K) = [ 

Jo 



T /r 1 1 

d s Qs + 2 Vp. ~ -a(p s )VH s 



VH s Jds 

• / T (d s K s G s - C{p s ){e G ° - 1) - A(p a )(e- d ' -l))ds. 



o 



for any (G, H) smooth functions. 

For G and H given by (15. 2p . we are going to check that 

(5.5) Jg,h(p, Q, K) ^ sup Jq h(p, Q, K). 

G,H 

Indeed, take H = H + F and any G, then it is easy to see that 

J &j6 (p,Q,K) = J 6yH { Pl Q,K) - ~y ds (a(p s ) |VF S | 2 > ^ J^ H (p,Q,K). 



18 THIERRY BODINEAU AND MAXIME LAGOUGE 

Moreover, if G = G + F and H any regular function, then from the expression of Jq (12. 9 j) 
and the identity ( 15.2ft we get 

J d6 (p,Q,K) = J G6 {p,Q,K) + / ds (C(p s )e G °(l + F s -e F °)) 

Jo 

+ [ ds(A(p s )e- G '{l-F s -e-^))^J GA (p,Q,K), 
Jo 

where we used that exp(a;) — x — 1 > for any igR. This completes the proof of (15.5ft . 
Finally for G and H given by (15.2ft . one can rewrite (15.4ft 

Ja#(p,Q,K) = l -f* (^pj (d s Q s + ~V Ps ) 2 Jds + £ (C(p s ) (l-e G ° + G s e G °))ds 

+ / (A(p s ) (l-e^ -G a e~ G ')) ds. 

This leads to the explicit form of the functional ( 12.13ft . □ 

Remark 5.2. For G and H as in ( 15. 2\) , then the conservation equation dtp = dtK — VdtQ 
implies that p obeys the hydrodynamic limit (13. 2p of the modified dynamics ^q H (13. ip . 

We introduce now the set of regular trajectories 

Definition 5.3. Denote by S the set of trajectories (p, Q, K) satisfying I(p, Q, K) < oo 
and such that 

• p is bounded away from and 1: there is e > such that e < p(t, x) < 1 — e for 
any (t,x) in [0,T] x [-1, 1]. 

• There exists two regular functions G and H such that (p, Q, K) is a weak solution 
(in the sense of (17. ip ) of 

(5.6) 

d t p(t,x) = ~Ap(t,x) - V (a(p(t,x))VH(t,x)) + C(p(t, x))e G ^ - A(p(t, x))e~ G ^ x \ 
p(t,±l)=p ± , pit = 0,x) = 7(2) , 

• For any smooth test function ip in C 1,1 ([0,T] x [—1, 1]) 

(5.7) +i jf (p s V(f s )ds+ jf (a(p s )VH s ^ s )ds, 

(5.8) (K T <p T ) 



£ (k s ^s\ ds = £ ds < [C(p s )e G ° - A(p s )e- G °] <p s ) . 
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The same argument as in Lemma [5.11 implies that for any trajectory (p, Q, K) in S 
I (p,Q,K) = J G>H (p,Q,K) 



(5g) = l -j* (a{p s )\VH s \ 2 )ds + £(C(p a ) {l-e G °+G s e G *))ds 

+ f (A(p s ) (l-e" Gs -G s e- G °))ds. 
Jo 

We prove now the lower bound for trajectories in S. Theorem 12 .41 is a direct consequence 
of relation (I5.9P and of the following Proposition. 

Proposition 5.4. For all open set O e £, 

liminf ilogP" W,qr,K») 6 O] >- M jg ons I( P ,Q,K). 

Proof. Let (p,Q,K) be in Of]S and satisfying I(p,Q,K) < oo. There is regular (G,H) 
for which p is a weak solution of (15.61) . Thanks to Lemma \3. 11 there is e > such that the 
trajectories concentrate in the set Bs tE (14.31) 

(5.10) limsup^logP^^J < -1. 

The function J G ,n,e defined in ( 14. 4p is continuous on £. Moreover, since po is bounded 
away from and 1, the function <p — > h ljPo (ip(0, •)) is continuous on D([0, T], Alo)- Let 
V C be an open neighborhood of (p, Q, K) such that 

V(<p,U,L)eV, \J G , Hte (<p,U,L)-J G ,H,e(p,Q,K)\<6, 

\h^, P0 ((f ) - h^ Po (p )\ = \h l!P0 (ip ) - h^(p )\ < 5. 
Using the change of measure (14.11) 

/ dF N 

^ exp{-iV(J Gi H,e(p,g,^) + /i 7 (po) + G ,H(iV- 1 , e ,5))}P^ GiH (Vn J B (5i£ ) . 
Thanks to (15. 10ft and the hydrodynamical limit for the perturbed process (13. 2\\ 

lim F% iG>H (V n B s , e ) = 1 , 

where we used the uniqueness of the weak solution of ( 15.61) (see the Appendix) to conclude 
that the probability of V converges to 1. This leads to 

liminf 1 logP^(O) > -J G;H , e (p, Q, K) - h 7>po (p ) + G)H (5, e). 

Letting e J, and <5 | 0, and since (G, if) satisfy ( 15.61) . the Proposition is completed. □ 
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5.2. Approximation for general trajectories. To complete Theorem 12.61 for general 
open sets, it remains to prove that 

Lemma 5.5. We assume (LI, L2). For any (p,Q,K) such that I(p,Q,K) < oo there is 
a sequence (p^ n \ Q^ n \ K^ n >) in S converging weakly to (p,Q,K) such that 

(5.11) I(p, Q, K) = lim I(j>W,QW, K™) . 

n— yoo 

Lemma 15.51 implies that 

(5.12) inf I(p,Q,K)= inf I{p,Q,K). 

Combining this identity and Proposition 15.41 proves Theorem 12.61 

5.2.1. Bounding the density away from and 1. We first approximate the density by tra- 
jectories bounded away from and 1. 

Lemma 5.6. Let P = (p, Q, K) be a path such that I(P) < +oo. There is Ps = (ps, Qs, Ks) 
with density ps uniformly bounded away from and 1 which converges to P and such that 

I{P) = \imI{P 5 ). 

<5->0 

Proof. Using Assumption (LI), we first establish a property of the functional Jo- We use 
the decomposition (15. ip of I . I\ is a convex functional of (p, Q) as it is the supremum of 
Jjj which are convex functional of (p, Q) (we used that a(p) = p(l — p) is concave). J 2 is 
not convex, but we use a trick introduced in [JLVj and decompose I2 as 

I 2 (p,K)= sup J 2 G (P,K)+ ( T dt(C(p t )+A(p t )). 
Gee 1 ' Jo 

Since A and C are concave, sup G J@(p, K) is a convex functional of (p, K). Thus the large 
deviation functional can be decomposed into two terms 

(5.13) J (p, Q, K) = J (p, Q,K)+ I dt (C{ Pt ) + A{p t )) , 

^0 

and Iq is a convex functional of (p, Q, K). We deduce that I is lower semi-continuous for 
the weak topology. 

We turn now to the approximation procedure. Let P = (<p, U, M) be the solution of the 
hydrodynamic equation (12.61) 

d t <p(t,x) = ±A(p(t,x) + C((p(t,x))-A((p(t,x)) 
d t U(t,x) = -§V<p(t,x) 
d t M(t,x) = C(<p(t,x))-A(<p(t,x)) 

with boundary conditions 

y?o(0, x) = 7 (x) , U{0, x) = 0, M(0, x) = 0, <p(t, -1) = p_, <p{t, 1) = p+ . 
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By construction I(P) = 0. We set P$ — (1 — 8)P + SP which has a density bounded away 
from and 1. As Io(P) < oo, the convexity implies that Io(Ps) ^ (1 — 5)Io(P) + SIo(P) 
so that 

limsup/o(A) <Iq(P). 

As Pg weakly converges to P, the lower semi-continuity of Iq implies 

J (P) ^liminf Io(P s ). 

<5->0 

Thus Iq(Ps) converges to h(P)- Finally /i 7 and p — > J"„ T dt (C(p t ) + A(p t )} are continuous 
for || • | |oo. This completes the Lemma. □ 

5.2.2. Time regularisation. We will prove 

Lemma 5.7. For any path P = (p, Q, K) such that I(P) < oo, there is a sequence regular 
in time P £ = (p t , Q e , K t ) converging weakly to (p, Q, K) such that I{P e ) converges to I{P)- 



In the following, only the regularity of K is needed in order to construct a drift G 
adapted to the non-conservative current (15.151) . However, the regularizing sequence has to 
satisfy the conservation law (12.101) so that p,Q,K will be approximated simultaneously. 

Proof. The proof is based on a time convolution and follows similar steps as in Lemma [5T6l 
We just recall the salient features of the proof (see [KL] for further details). 

Let ij) e be a C°° approximation of unity such that ip e = outside [0, e] and Jjf i/; e = 1. To 
take the convolution product of (p, Q, K) with ip e , we have to extend the path (p, Q, K) 
beyond the time T. Let r be the solution of the hydrodynamic equation (I2.6P 

d s r s = ^Ar s + C(r s ) - A(r s ), with r = Pt ■ 



2 



We set for s > T 



p s = r s _ T , Q s = Qt ~ \ f du Vp u , K s = K T + [ du \C(p u ) - A(p u )} . 

£ Jt Jt 

Now we define 

{p e (t,x) = p*i/j £ (t,x), 
Q e (t, x) = Q * i/> e (t, x) - Q * V> e (0, x), 
K e (t, x) = K * ifj e (t, x)-K* ^ e (0, x) . 

The path (p e ,Q e , K e ) satisfies the relation (12.101) and has initial currents equal to 0. We 
use the decomposition (15.11) and approximate independently the functionals 1% and I2 by 
using convexity properties deduced from Assumption (LI). □ 
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5.2.3. Non regular drifts. Thanks to Lemmas 15.61 and 15.71 it is enough to consider a tra- 
jectory (p, Q, K) regular in time with a density uniformly bounded away from and 1. We 
are going to associate to (p, Q, K) the drifts H, G as for the trajectories in S introduced 
in Definition 15.31 Note that the drifts H, G can be non regular in space. 

As d s K exists, the drift G is defined (as for the regular trajectories) as the solution of 

(5.15) Va: G [-1, 1], d s K(s, x) = C(p(s, x))e G{s ' x) - A(p(s, x))e~ G{s ' x \ 
Note that K, G solve the equation (15. 8p and as in the regular case (15.91) 

(5.16) h{p,K)= f (C(p s ) (l-e G °+G s e G °)) + (A( Ps ) (l - e~ G ' - G a e~ G ') ) ds. 

Jo 

The functional I\ is the same as the functional of the SSEP (without Glauber rates). 
Thus the drift of the conservative dynamics can be approximated thanks to the Riesz 
representation theorem as in [KLl IBDGJL2] . For (p, Q) such that Ii(p,Q) < oo, there is 
H in L 2 ([0,T],H!(] - 1, 1[)) such that 

(5.17) h{p,Q)= f T r a(p( ^ x)) (VH(s,x)) 2 dxds. 

Jo J-i 2 

Moreover (15. 7p holds. Note that the time regularity is not needed to derive (15.171) . 

Finally, we check that p is a weak solution of (17. ip with the drifts G, H . As the large 
deviation functional is finite, p belongs to A so that Q(p) is finite and (17. 2p is satisfied. 
Combining the conservation law (I2.10p with (15. 7p and (15. 8p shows that (17. 3p holds. 



5.2.4. Approximation paths in S. We approximate now H, G by regularized drifts H n ,G n 
and prove that the associated sequence of paths (p( n \ KM) in S converges to (p, Q,K). 

Let H n be a sequence of C°° functions converging to H in L 2 ([0, T], Hi(] — 1, 1[)). Since 
C and A satisfy assumption (LI) and p is uniformly bounded away from and 1 (Lemma 
15. 6p . C and A are either uniformly bounded away from 0, or uniformly equal to 0. We 
will consider only the case A > 0, C > as the other case follows in the same way. From 
(15TTB]) . we get that |G|e |G| e ^([O.T] x [-1, 1]), so that there is G n such that 

(5.18) / / \G n e ±Gn - Ge ±G \ >0 and / / \e ±Gn - e ±G \ > 0. 

Jo J-i n ~ > °° Jo J-i n ^°° 

Let (p( n \ Q^ n \ K^) be the sequence of paths in S associated to the regular drifts H n , G n . 
In particular the weak solution of 

( . ! ()) ■ dtpW = ^ApW - V(a(p( n ))V#„) + C(p( n ))e G ™ - A(p( n ))e" G " 

p(")(t,±l 2 )=p ± , pW(0,-) = p(0,-), 
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and Q (n) and are defined as 

(5.20) Q {n \t,x) = f (-l-Vp {n \s,x) + a(p {n) (s,x))VH n (s,x))ds, 

Jo * 

K (n \t,x) = [\c(p (n \s,x))e G " {s ' x) -A{p {n \s,x))e- G " {s ' x) )ds. 
Jo 

From energy estimates (using similar bounds as in the Appendix (17.81) ) we get 

(5.21) J dtj dx{Vp {n) f < C, 

where C is a constant independent of n. As p^ n ' and V ' p^ are bounded sequences in L 2 , 
they are tight in the weak topology. We want to check the uniqueness of the limiting 
points. Let R be a limit point of p( n \ We will prove that R is a weak solution of (17.11) . 
From Section I5.2.3[ we know that p is also a weak solution of (17.11) so that the uniqueness 
of the weak solutions (see the Appendix) will imply that p = R. 

We want to take the limit in the weak formulation of (I5.19p . By construction, for any 
smooth function (p on [—1, 1], we have 

ds I dxp {n) {s,x)Vp{x) = \p+<p(l) -p-(p(-l)] - [ ds [ dxVp {n) {s,x)p{x). 



-i 



By weak convergence of the subsequence, one has 

ds dx R(s, x)V(f(x) = \p + (p(l) — p_(p(— 1)] — ds dxVR(s,x)(p(x). 
Jo J-i Jo J-l 

Furthermore (15.21ft implies that the limit V-R is also in L 2 . It remains to take the limit in 
n in the equation 

(5.22) T p(™)(% = / IV ~\f [ Vp^Vif + f [ Vipa{p^)VH n 
J -i J-i 1 Jo J -i Jo J -i 

+ f fc(p (n) )e G ^- f I A(p (n) )e- G 'V- 
Jo J-i Jo J-i 

The difficulty is to treat the non linear terms. We proceed term by term and start with 
the non-linearity C in ( 15.221) (the term in A can be controlled in the same way) . By (15.181) , 
exp(G n ) converges to exp(G) in L 1 thus it is enough to check that 



(5.23) lim f I C(p {n) )e G p = f C C{R)e G ^> 

n ^°°Jo J-i Jo J -i 



24 



THIERRY BODINEAU AND MAXIME LAGOUGE 



We write 



(5.24) f f {C(p^)-C{R))ipe G = f f (<7(p<»>) - C{p™ * n))<p e G 
Jo J-i Jo J-i 

+ f I (C(p {n) *L S )-C(R*Ls)) V e G + f I (C(R*i 5 )-C(R)) V e G . 
Jo J -i Jo J-i 

As 5 goes to 0, R * l§ converges to R. The convergence in the weak topology implies that 
p^ * is converges a.s. to R * l$ when n goes to infinity. Suppose that 



(5.25) 



lim sup 



t rl 



J-l 



C(p^)-C(p^*t s ))e^ 



(n) 







then choosing 5 small and then n large, the convergence in (15.231) follows. 
To derive (I5.25p . we first note that for x G] — 1 + 5, 1 — S[ 



(5.26) p (n) * i s (x) - p {n \x) 

dyL S (x-y)[pW(y)-pW(x)} 



-i 



dy l s (x - y) I du d u p {n) (u) 



du(d uP W(u)) 2 / dvL 5 (x-yW\x-y\^2VsJ / du{d uP ^{u)f. 



Using the uniform bound (15.211) there is a constant C independent of n such that 



(5.27) 



7 

J-i 



(p (n) * L 5 - p {n) ) 2 < CS. 



It remains to check (15.251) . Let m be the uniform measure on [0, t] x [—1,1] given by 
^dsdx. Given ip(s,x) a non-negative function on [0,t] x [—1, 1] such that m(^) = 1, then 
the entropy inequality implies that for any function ip 

(5.28) m(<p(s, x)ip(s, x)) < m(ip(s, x) log (ip(s, x))) + logm(exp (<p(s, x))) . 
Since A and C are positive then m(e' G '|(jr|) is bounded and (15. 28ft implies 



m 



,|G| 



me 



^ Esm 



,\G\ 



|C7(pW*^)-C(pW)|M 

e |G| 

log r- 

v m(el G l) m(el G l) 

where e$ vanishes (in a suitable way to be determined later) as S goes to 0. Recall that 



e s \ogm exp (-|C(P (n) * is) - ^(P W )IM) 
£5 



VxG[0,M], exp(x) < 1 + x + exp(M)x 2 and log(x) < x - 1 . 
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As \C(p^ * ig) — C (pW)\\(p\ < 2 1| C || oo || V 3 1| oo 5 the previous inequalities imply 
logm fexp (-\C(pM * is ) - C(pW)\\tp\) 



^5 



1 /H|C||oo|M 



^ 3 exp 



m (\C(pW * Ls ) - C(p (n) )IM + (|C(pW * l s ) - C7(pW)||^|) : 



£ 5 \ £ S 

We choose es such that 

r n r /l|C||oo||^||oo\ n 
lim £5 = (J, hm exp = . 

5^0 eg \ £5 J 

Combining the previous inequality and (I5.27p . we conclude that f!5.25j) holds. 
We turn now to the non-linearity a. Using the decomposition 

f I [o{p {n) )-a{R)]V^ s VH s = f f [<t( p W) -<r(pW*L S )]V<p s VH 8 
Jo J-l Jo J-l 

+ 11 [°(P {n) *is)-cr(R*Ls)]V^ s VH s + [ [ [a(R * t s ) - cx(R)} V<p s VH S 
Jo J-i Jo J-l 

the last two terms vanish in the limit as in ( I5.24p . The first term can be controlled thanks 
to the uniform bound (15.271) and a Cauchy-Schwartz estimate as VipVH belongs to L 2 . 
This concludes the convergence of (I5.22p . 

Following the same proof, (15T20D implies that (Q^ n \ K^) converges weakly to (Q,K). 

5.2.5. I— convergence. We finally complete the proof of Lemma [5.51 Let (p^ n \Q^ n \ K^) 
be the regularizing sequence defined in ( I5.19p . We start by proving 

(5.29) limsup/o(p (ri) ,g (n) ,K (n) ) < I {p,Q,K). 

n—>oo 

From the expression (I5.17P of I\ and the weak convergence of p^ to p, we get 
I x {p, Q) = C C a{p ^ x)) VH(s,x fdxds = lim C C *#VF( S ,)ta 



jo J-i 2 s ^°Jo J-i 2 

= lim lim - / / V//(s,x) 2 ( x([p( n) * l s ](x)) dxds 

^ lim lim sup ~ [ [ (yH(s,x)) 2 [o-(p^)*i s (x)]dxds, 

n^oo ^ Jo J — 1 

where we used the concavity of o in the last inequality. Notice that 

»T i-l j-T rl 

(\7H(s,x)) 2 [o-(p^)* L 5 {x)]dxds = / / (Vif(s, x)fa{pf\x)) dxds 

u J-l ' " Jo J-l 

T f 1 r 1 

/ (VH s ) 2 *i s (x) - (VH(s,x)) 2 a{pf\x))dxds. 
J-i L -I 
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As a is bounded and (VH) 2 is integrable, the last term vanishes uniformly in n when 5 
goes to 0, so that 

h{p,Q) > limsup^ I f VH(s,x) 2 a(p^\x))dx ds. 

As converges to \7H in L 2 and a is bounded, we conclude that 

(5.30) h(p,Q) > Mrnsnp h(p (n \Q (n) ) . 

Similarly for I2, the concavity and boundedness of C and A imply that 
I 2 {p, K) > lim sup f [ C(p (n) ) (1 - e G + Ge G ) + A{p {n) ) (l - e~ G - Ge" ) . 

Then we use the convergence properties on G^ n > (I5.18P and the fact that C and A are 
bounded to conclude that 

(5.31) Hp,K) ^ limsup/ 2 (p (n) ,is: (n) ) 

n— >oo 

Estimates f l5T30|) and f lOTj) imply (15^9|) . 
The converse inequality 

(5.32) hminf 7 (p (n) , Q (n) , ^ (n) ) >I (p,Q,K), 

n— »oo 

can be deduced from the weak convergence of (p (n) ,Q (n) , if (n) ) to (p,Q,K) and from the 
decomposition ( 15.1 3ft : the limit follows for the term I from its lower semi-continuity and 
the convergence of the second term / Q T dt (C(p5 n) ) + A(p| n) )) can be obtained as in (15.231) . 

6. The density large deviations 

In this section, we recover the density large deviation principle (first derived in [JLVj ) 
by optimizing the functional I over the currents. We assume that the rates A, C sat- 
isfy assumptions (LI, L2). The contraction principle |DZj implies that the density large 
deviation functional is given by 

(6.1) F{p) = inf Up, Q, K) , 

where the infimum is taken over the currents (Q,K). Using the approximation procedure 
of section I5.2[ we will check that it is enough to consider regular density profiles and the 
modified functional 

(6.2) f(p)= inf I(p,Q,K), 

(Q,K) teg 



where the infimum is taken now over the regular currents (Q, K). 
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Step 1. In this first step, the explicit solution of the variational problem ( 16. 2 p is computed 
for regular trajectories. The functional (16. 2p can be rewritten as (see (15.91) ) 

P{p)= wt{\ f \a{p s )\VH s \ 2 )ds + f T (C(p s ) (l- e G ° + G s e G °)) ds 



(6.3) 



+ J\a{ Ps ) (l- e - G °-G s e- G °)) dsj 



where G, H are smooth functions such that the conservation relation ( I2.10p holds 
(6.4) 8 sPs = ^A Ps - d x (a(p s )VH s ^ + {C(p s )e G ° - A(p s )e- G °). 

We first check that the infimum is reached for functions H with boundary conditions 
H(s, —1) = H(s, 1) = for any time s. Let / be a smooth function in [0, T\. Perturbing 
H into H(s, x) + f(s) J* -^-h u ^ du, we see that the conservation law (16. 4p is preserved (it 
simply amounts to adding a constant conservative current) and 



T I pi \ 2 



a(p(s : x)) 



^ds = jf (trips) \VH s \ 2 )ds 

As H minimizes the integral, this implies that H vanishes at the boundaries. 

Suppose that an extremum is reached at (if, G) and consider a perturbation with the 
new drifts H + h and G + g. Then the constraint (16. 4ft implies that h, g satisfy the relation 

(6.5) <9 x (a( Ps )V/i s ) +g s (C(p s )e G ° + A{p s )e~ G °) = 0. 
A perturbation of (16.31) around the extremum H, G leads to 

(6.6) 0=/ (a{p s )VH s Vh s ) ds + I (g s G s (C(p s )e G ° + A(p s ) e - G °)) ds . 

Jo Jo 

Since H vanishes at the boundaries, the relation ( 16.51) combined to ( 16.61) leads to 

= f (g s [G s - H a ] (C(p s )e G ° + A{p s )e- G °)) ds . 
Jo 

This holds for any g so that the extremum is such that G = H with H determined by 

(6.7) d tP = X -Ap - V (a(p)VH) + C(p)e H - A(p) e - H . 

Thus if H satisfies (16. 7p and G = H, then an extremum is reached for the corresponding 
currents (Q,K). Since the functional I(p,Q,K) is convex wrt (Q,K) (thanks to the 
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representation (12.81) . (I2.9p ). the extremum (Q,K) has to be a global minimum, thus 

Hp) = \ [ T ( ( r(p s )\VH s \ 2 )ds+ f (C(p s ) (l-e H ° + H s e H °)) ds 
* Jo Jo 

+ f (A(p s ) (1 - e~ H * - H s e- H °)) ds . 
Jo 

Step 2. To approximate J 7 (16. ip in terms of T (16. 2p . we first check that the minimum 
is reached in (16.11) . Consider a sequence (Q n ,K n ) which realizes the infimum. By the 
tightness argument (section 14.21) the sequence belongs to a compact set and therefore has a 
weak limit (Q*, K*). From the lower semi- continuity of the functional (Q, K) — > I(p, Q, K), 
this weak limit is a minimizer 

(6.8) F{p) = I{p,Q\K*). 

Let (p e , Q*, K*) be a regular sequence (as in (I5.19P ) converging to (p, Q*, K*) such that 

I{p e ,Q*,Kt)^I{p,Q*,K*). 
As P{p £ ) < I(p e ,Q*,K*), one has 

limsup^(p £ ) < \imI(p ei Q* £) K*) = JF{p) . 

e^O 

Using (16.71) . there are regular (Q e ,K e ) such that 

f(p E ) = i(p £ ,Q e ,k e ). 

The sequence (p e , Q e , K £ ) has a bounded large deviation cost and thus it belongs to a 
compact set. There is a subsequence such that (Q £ ,K £ ) converges weakly to (Q,K). One 
gets 

\imMf(p £ )>I(p,Q,K)>F(p). 

e^O 

This limit follows from the decomposition (15.131) . The term J converges by lower semi- 
continuity and the convergence of the second term J Q T dt(C(p £ (t)) + A(p £ (t))} can be 
obtained as in (15.231) (p £ is just a subsequence extracted from p^). 

Combining both estimates, we deduce that T can be approximated by JF. 

7. Appendix : Uniqueness of the weak solutions 

Let H be in L 2 ([0, T], Hi(] - 1,1[)) and \G\ exp(|G|) in L x ([0,T]x] - 1,1[). Given an 
initial data 7, a weak solution of 

(7 ^ f d t p = - d x (a(p)d x H) + C(p)e G - A(p)e~ G 

K } \ P(t,±l) = p±, p(0,x)= 1 (x) 

is defined as : 
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• The density p is in L 2 ([0,T], H^Q - 1, 1[)), i.e. there is a function in L 2 ([0,T]x] - 
1,1[) which will be denoted by Vp such for every t G [0, T] and every function 

(7.2) ds dxp(s, x)V<p{x) — {p+<p(l) — p_<p( — l)}t = ds dxV p{s,x)f{x) , 
Jo J-i Jo J-i 

where p± are fixed boundary conditions. 

• For every t G [0, T] and every function <p G C 1 ([— 1, 1]) vanishing at the boundaries, 

(7.3) / dxp(t, x)ip(x) — / dx / y(x)(p(x) = — ds dxVp(s,x)Vip(x) 



t ,i 

ds I dx u(p(s, x)) VH(s, x) V<p(c 
+ I ds I dx(C(p(s,x))e G ^ s ' x) - A(p(s,x))e- G{s ' x) ) <p(x) . 



ft /.l 

(is / 

'o i-i 

The hydro dynamic limit (12.61) corresponds to H = G = 0. 

In this Appendix, we derive the uniqueness of the weak solutions. The main technical 
difficulty comes from the fact that G,d x H are unbounded (see [E] for bounded drifts). 
Note that at this stage Assumption (LI) is irrelevant. We will rely on Assumption (L2) 
on A, C which can be interpreted as follows. Equation (17.11) in a strong form reads 

d t p(t, x) = Ap(t, x) - d x (a(p(t, x))VH(t, x)) - V^ x) (p(t, x)) , 

where the reaction term is determined by the space-time dependent potential Vc(t,x) with 
Vg(p) = —e 9 C(p) +e~ g A(p). Assumption (L2) ensures that the potential Va(t,x) is convex 
uniformly in G(t, x). Thus the reaction and the diffusion terms are both contractions and 
the solution will be unique. We adapt to our framework the argument of |LMSt IFLMj . 

We consider two initial datas pi, pi and the corresponding weak solutions p\p 2 . We 
are going to prove that the L 1 -norm \\p\ — p 2 ||i decreases in time. For a given 5 > 0, we 
introduce the regularized absolute value 

(7-4) U s {v) = ^1{M<<5} + (M - 2) k\v\>s} ■ 

Define 

(7.5) Vs = |(s, x) G [0, T] x [—1, 1], such that \p\s, x) - p 2 {s, x)\ < 5 } . 
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Step 1. We are going to check that for times t <t' 



(7.6) J dxU s (p l (t' » - p 2 (t' - J dxU s (y-(t,x) - p 2 (t,x)^j 

= -~ J J dsdxl^^v^Vip 1 - p 2 )(s,x) |v(p x - p 2 )(s,x) + a(s, x)VH(s, x)| 



+ J dsdxU' 5 (p 1 (s,x)-p 2 (s,x)){c(s 1 x)e G{s ' x) -A(s,x)e- G{s ' x) 
where we set 

A(t,x) = A(p l (t,x))-A(p 2 (t,x)), C(t,x) = C(p 1 (t,x))-C(p 2 (t,x)), 

a(t,x) = a(p l (t, x)) — a(p 2 (t, x)) . 

To prove (17.61) . we follow the regularization scheme introduced in [BLMJ (see the proof 
of their Theorem 4.6). For e > 0, denote by : [—1, l] 2 — > R + (resp R^) the resolvent 
of the Dirichlet Laplacian A_d (resp Neumann A^r) 

Rf = (Id - sA D )-\ R» = (Id - eAk)- 1 . 

The mollified trajectory is defined by 

p £ (t,x) = p(x) + R? (p t - p)(x) , 

where p stands for the linear profile between p + and p_. Note that the resolvent of the 
Dirichlet Laplacian preserves the boundary conditions of the mollified trajectory. 
As p is a weak solution of (17.31) . one has 

d t p £ (t,x) = d t J dyR?{x,y)(p(t,y)-p(y)) 

^ dy d y R°(x, y){-d y (p(t, y) - p(y)) + a(p(t, y))VH(t, y) } 

dyR?(x,y) (C(p(t,y))e G ^-A(p(t,y))e- G ^) , 



where we used the fact that R^(x, 1) = R^(x, —1) = for any x. 
From the relation d y R^(x,y) = —d x R^ (x,y), one has 

d t p £ (t,x) = d x J dyR?{x,y){d y (p{t,y)-p{y))-a(p{t,y))VH{t,y)} 
(7.7) + / dyR?(x,y) (C(p(t,y))e G ^ - A(p(t,y))e~ G ^) 



Taking the time derivative we obtain 



dtJ dxU s (p 1 > £ (t,x)-p 2 > e (t,x)) = j dxU' s (p 1 > e (t,x)-p 2 > £ (t,x)}d t (p 1 > £ (t,x)-p 2 > e (t,x)) 
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As d x R?(x,y) = -d v R*(x,y), we get 



1 



d x (p 1 > £ {t,x)-p 2 <%t,x)) = - I ^dyd y R?(x,y)(p\t,y)-p 2 (t,y)) 

dyR?(x,y)(d y p 1 (t,y)-d yP 2 (t,y)). 

where we used that p 1 , p 2 are in L 1 ([0, T], Hi(] — 1, 1[)). Thus we can write 

d x [U' s (p 1 '%t,x)-p 2 > £ (t,x))] = l {x& vi t }Re (d x {p x -p 2 )){t,x) , 
with the notation 

Vgj = (x G [— 1, 1], such that |p 1,e (t, x) - p 2 ' e (t, x)\ < 5 } . 

Using relation (17.71) . one obtains 

d t J dxU s (p 1 ' £ (t,x)-p 2 ' £ (t,x)') 

dx J dy l {x€V e t} R* [d x (p l - p 2 ))(t, x) 
R?(x, y) {d y (p 1 (t, y) - p 2 (t, y)) + a(t, y)VH(t, y) } 
+ // <kdyl7Z((p*-^)(f,s)^ . 



-l 

Taking the limit as e tends to 0, we recover (17.61) . 
Step 2. 

First note that from assumption (L2), one has 

^((p 1 -p 2 )(t,x)) (C(t,y)e G ^-A(t,y)e- G ^) < 0. 
Thus the reaction term acts as a contraction. As a(p) = p(l — p), then 

V(s,x) G Vs, \a(s,x)\ < 5. 

Thus (17.61) implies 



t' 



(7.8) (U s {pl, - p 2 )) - (U s (pl - p 2 )) < -- J ds(l {ix>s)eVs} (Vpl - p 2 ) > 

+ J* ds(l {{X!s)eVs} \V(pl- p 2 s )\\VH s \) , 



using the fact that 



2\V(p 1 - p 2 )(s,x)\ \VH{s,x)\ < -IVip 1 -p 2 )(s, x) \ +8\VH(s, x) 

o 
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we obtain 



(U 5 (p l t ,-p 2 t ))-{U s (p\-pl)) < ~jT ds(l {[xAeVs} (Vpl-p 2 s ) 2 ) + 6 J* ds(\VH s \ 2 

<5^ ds(\VH s \ 2 ). 

Recall that H belongs to L 2 ([0, T], Hi(] - 1, 1[) thus as 5 tends to 0, the LHS converges to 
0. Furthermore Ug converges to the absolute value function. This implies that 



Vt < t' 



dx 



-i 



p L (t',x)- p 2 (t',x) 



< 



dx 



-i 



p\t,x) - p 2 (t,x) 



from which the uniqueness of the weak solutions follows. 



□ 
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